Quantum effects with an X-ray free electron laser 
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A quantum kinetic equation coupled with Maxwell's equation is used to estimate the laser power 
required at an XFEL facility to expose intrinsically quantum effects in the process of QED vacuum 
decay via spontaneous pair production. A 9TW-peak XFEL laser with photon energy 8.3 keV 
could be sufficient to initiate particle accumulation and the consequent formation of a plasma of 
spontaneously produced pairs. The evolution of the particle number in the plasma will exhibit non- 
Markovian aspects of the strong-field pair production process and the plasma's internal currents 
will generate an electric field whose interference with that of the laser leads to plasma oscillations. 

PACS numbers: 42.55.Vc, 41.60.Cr, 12.20.-m, ll.15.Tk 



X-ray free electron laser (XFEL) facilities are planned 
at SLAC Q: namely the Linac Coherent Light Source 
(LCLS), and as part of the e+ linear collider project 
(TESLA) at DESY [||. They propose to provide narrow 
bandwidth, high power, short-length laser X-ray pulses, 
with good spatial coherence and tunable energy. It is an- 
ticipated that the realisable values of these parameters 
will enable studies of completely new fields in X-ray sci- 
ence, with applications in atomic and molecular physics, 
plasma physics, and many other fields |Q . 

A unique ability of these facilities is to provide very 
high peak power densities. For example, a, P — 0.2 TW- 
peak laser at a wavelength of A = 0.4 nm, values which 
are reckoned achievable with current technology ||^ , can 
conceivably produce a peak electric field strength 

Ea = VmcP/iirX^) « 1.2 X IQi*^ V/m. (1) 

Boosting P to 1 TW and reducing A to 0.1 nm, which 
is theoretically possible [^, would yield an order-of- 
magnitude increase: Eg = 1.1 X 10"V/m. Electric 
fields of this strength are sufficient for an experimental 
verification of the spontaneous decay of the QED vac- 
uum |,|||10. 

It is a long standing prediction that the QED vacuum 
is unstable in the presence of a strong, constant elec- 
tric field, decaying via the production of e+ pairs 
In such fields, appreciable particle production is certain 
if the strength exceeds i?cr := ™e/s — 1-3 x 10^*V/m. 
(We subsequently use h ~ c ~ 1.) The proposed XFEL 
facilities could generate E w 0.1 -Ecr- (NB. Here "con- 
stant" means that the field must be uniform over time- 
and length-scales much greater than the electron's Comp- 
ton wavelength: l/rUe ~ 0.4 pm.) 

A single laser beam cannot produce pairs (For a 
light-like field F^^F'^'^ — and hence the vacuum sur- 
vival probability is equal to one.) Nevertheless, if two 
or more coherent beams are crossed and form a standing 
wave at their intersection, one can hypothetically pro- 
duce a region in which there is a strong electric field 



but no magnetic field. The radius of this spot volume 
is diffraction limited to be larger than the laser beams' 
wavelength: > A, and the interior electric field could 
be approximately constant on length-scales approaching 
this magnitude. The period of the electric field is also de- 
termined by A. Hence at an XFEL facility one might sat- 
isfy the length-scale uniformity conditions noted above. 

However, the planned laser pulse duration: Tcoh ~ 
100 fs, is large compared to the laser period: ~ 1 as, 
and thus the time-dependence of the electric field in 
the standing wave may materially affect the pattern 
of observed pair production; i.e., vacuum decay might 
be expressed via time-dependent pair production. This 
possibility can only be explored using methods of non- 
equilibrium quantum field theory pO| or an equivalent 
quantum kinetic theory ijll] , [l^ , ^ . The latter procedure 
was employed in Ref. j7| , wherein it was shown that pair 
production occurs in cycles that proceed in tune with the 
laser frequency. While that does not lead to significant 
particle accumulation, the peak density of produced pairs 
is frequency independent, with the consequence that sev- 
eral hundred pairs could be produced per laser period. 

The proposed XFEL facilities offer a first real chance 
of observing the decay of QED's vacuum, a profound and 
nonperturbative quantum field theoretical effect. Never- 
theless, with a quantum Vlasov equation, one can ask 
for more. This equation yields the time-dependence 
of the single particle distribution function: f{p,t) := 
(at(t) ajr(t)), and hence can be used to estimate the laser 
power required to achieve an accumulation of e~ e"*" pairs 
via vacuum decay. Furthermore, in quantum field theory 
the particle production process is necessarily non-local in 
time; i.e., non-Markovian, and dependent on the parti- 
cles' statistics. These features are preserved by the source 
term in the Vlasov equation. (The Schwinger source term 
is recovered in a carefully controlled weak-field limit |p2[.) 
Consequently, one can identify the laser parameters nec- 
essary to expose negative energy elements in the particle 
wave packets. The system exhibits this core quantum 
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field theoretical feature when the time between produc- 
tion events is commensurate with the electron's Compton 
wavelength. Along with accumulation comes the possibil- 
ity of plasma oscillations, generated by feedback between 
the laser-produced electric field and the field associated 
with the production and motion of the e+ pairs; and 
also collisions. (These features are reviewed in Ref. [p^.) 

We model an "ideal experiment," and assume r„ — \ 
and a nonzero electric field constant throughout this vol- 
ume while the magnetic field vanishes identically. This is 
impossible to achieve in practice and therefore the field 
strengths actually achievable will be weaker than we sup- 
pose. Hence our estimates of the laser parameters will be 
lower bounds. Our model is represented by 
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^cxt(i) = (0,0, £^0 sinm), ^ 



:27r/A. (2) 



Particle production occurs throughout the coherence 
spike length, Tcoh- However, collisions and radiation may 
become important as time progresses and therefore, to 
simplify our analysis, we focus on the first 100 laser pe- 
riods; i.e., the first ^ 0.1% of the laser pulse. 

In spatially homogeneous fields the kinetic equation 
has the form 



dt 



= S{p,t) + Cip,t), 



(3) 



where S{p^ t) is a source term describing particle pro- 
duction and C{p,t) is a collision term, which properly 
includes annihilation. The form is intuitive: the rate of 
change in the single particle distribution function, /(p, i), 
is determined by competition between particle produc- 
tion, and collisions and annihilation. The source term 
can be calculated directly using quantum mean field the- 
ory pT| , p^ . The collision term, however, is more com- 
plicated and not directly accessible in the mean field 
approximation. It only becomes important, though, in 
dense plasmas and because we deliberately avoid that 
situation, we proceed by setting C = 0. 

The distribution of particles is described by f{p,t). 
These particles are accelerated by £^ext(i), and the associ- 
ated currents generate an opposing electric field, Eint{t)- 
It is the sum: E{t) — Ecxt{t) + Eint(,t), which is finally 
responsible for particle production. Thus, to properly 
represent the system's evolution, one must solve a Vlasov 
equation coupled with Maxwell's equation Q: 



dfip,t) 
dt 



Eint it) 



eE{t)el eE{t')[l^2f{p,t')] 
2u;^{p,t)Jo u;Hp,f) 
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(2^)3 L0{p,t) 

uj^ip,t)dfip,t) ^ eEjt) ej 
eE{t) dt 8uj^{p,t) 



(5) 
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FIG. 1: Number density calculated with Eq = 0.5 Ecr- Solid 
line: solution of Eqs. (^); dotted line: solution obtained 



using a low density approximation to S{p,t), Eq. 
oscillations are tied to the laser frequency, see Eq. 
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with: p = (p±,P3); ei = ml + p\; uj^{p,t) = t 
P||(t) =P3-eA(t). 

Quantum statistics affect the production rate through 
the term "[1 — 2 /]" in Eq. (^, which ensures that no mo- 
mentum state has more than one spin-up and one spin- 
down fermion. In addition, both this factor and the "cos" 
term introduce non-Markovian character to the system: 
the first couples in the history of the distribution func- 
tion's time evolution; the second, that of the field. 

There are two control parameters in Eq. (^): the 
laser field strength, Eq, and the wavelength, A. We fix 
A = 0.15 nm, which is achievable at the proposed XFELs. 
(NB. By assumption, the volume in which particles are 
produced increases with A^ whereas the field strength de- 
creases with 1/A: there is merit in optimising A.) Our 
study will expose additional phenomena that become ob- 
servable with increasing Eq- 

The particle number density: 



n{t) = 2 



d^p 
(2^ 



f{p,t). 



(6) 



is plotted in Fig. Q for E = 0.5 Ec - Our results for 
T ^ t/X ^ 100 are accurately fitted by 

n{T;Eo) - ao{Eo) sin^ 2TrT + p{T, Eg) T , (7) 
piT,Eo) - piEQ)+p'{Eo)T. (8) 

We can therefore use p(T,Eo) to quantify the rate of 
particle accumulation. This rate is very small for Eq ^ 
0.2 Ed- and, while noise in n{t) prevents a reliable nu- 
merical determination of p, p' , it is nevertheless clear that 
in this case lengthening Tcoh will not materially increase 
the number of particles produced. 

It is apparent from Fig. however, that the situation 
is very different for Eg — 0.5i?cr- The solid curve is 
described by Eq. (0), with uq — 1.2 x lO^^^fm"^^ p = 
5.4 X 10-12 fn^-3pgj.jQj-i p/^p ^ 0.0033/period. Clearly, 
the accumulation rate is approximately constant. 
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FIG. 2: Solid line: accumulation rate, p{0,Eo), Eq. M 
dashed line: Schwinger rate, Eq. (^); dot-dashed line: Eq. (| 
with as ~* 0.305, bs ^ b = 1.06. 



FIG. 3: Peak particle number density versus laser field 
strength. There is a striking qualitative change at Eo ~ 
0.25 i?cr, which marks the onset of particle accumulation. 



Solving the coupled system of Eqs. (|5|) numerically 
is straightforward but time consuming. We have repeated 
the procedure for a range of values of Eq and plotted the 
calculated values of p{Eo) in Fig. ^. These results can 
be used to obtain a rough estimate of the accumulated 
number density via: n(T, _Eo) ~ p{Eq)T. We also plot: 



PS = as 



Eo_ 

Ed 



^-bs.E^./Eo^ as 



1, b, 



1 : 



(9) 



i.e., the Schwinger rate. As can be anticipated from the 
energy budget, the constant-field Schwinger rate bounds 
our results from above. This comparison provides a con- 
text for the estimates in Ref. ||^. 

In Fig. ^ we plot the peak particle number density, 
n{t^), where = (4A: - 3)A/4 is the time at which E{t) 
is maximal during laser period no. k. The figure dis- 
plays a qualitative change in the rate of particle produc- 
tion with increasing laser field strength. [n(kX) w n(t^) 
if and only if there is significant particle accumulation, 
otherwise n{kX) <^ '^(^fc )•] 

That change is emphasised by Table |, which reports 



N{t>,Eo) ■.= anit>,Eo): 



(10) 



i.e., the peak particle number. For Eq < 0.25 Ea there 
is no significant accumulation of particles: N{tQ,Eo) « 
N(tfg,Eo), etc., just as observed in Ref. 0. However, 
for Eq > 0.25 Ecr there are more particle pairs after each 
successive laser period; e.g., Eq — 0.35 i?cr brings an or- 
der of magnitude increase in N over the first 100 laser 



TABLE I: Peak particle number, Eq. (1^), at selected values 
of the laser field strength, Eq in Eq. (2), with k = 50. 



Eq I Ec; 



0.1 



1.3 



0.2 



5.4 



0.3 



17.9 



0.4 



148 



0.5 



1110 



periods. At such values of Eq one is in a domain where 
numerous electrons and positrons produced in a single 
period are accelerated to relative longitudinal momenta 
that are sufficient to materially inhibit annihilation. This 
ensures that many pairs remain when the next burst of 
production occurs. Consequently 7V(i^,i?o) grows con- 
siderably with increasing k. (NB. Accumulation means 
collisions can become important and should be included 
in the kinetic equation if quantitatively reliable results 
are required. As we neglect these effects we do not follow 
the system's evolution into the high density regime.) 

In the low density limit: f(j),t) <^ 1, the distribution 
function disappears from the right hand side of Eq. (||), 
which then yields an algebraic solution: 



dT 



eE{T)e\ 



dt 



, eE{t') 



2u;^ip,t)Jo uj^{p,t') 



X cos 



2^ dr' cj(p,t') 



(11) 



We anticipated this in Fig. ||, which contrasts the low- 
density result with the complete solution for Eq = 
0.5 Ecr, which we now know generates significant par- 
ticle accumulation. The comparison emphasises that, 
with particle accumulation, quantum statistical effects 
in the source term are observable. These effects are es- 
sentially non-Markovian. The difference grows with in- 
creasing field strength because Pauli blocking becomes 
more important with increasing fermion density. 

We now return to Eq. (||); i.e.. Maxwell's equation: 
E ^ —j. The right hand side has two components: the 
first term, proportional to f{p,t), is a conduction cur- 
rent tied to the particles' motion; and the last two terms 
express a polarisation current that is tied to the pair pro- 
duction rate ||l^ . While these terms remain small the to- 
tal electric field is determined solely by the applied laser 
beams. However, in the domain of particle accumulation 
the internal currents can modify the electric field in the 
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FIG. 4: Ratio of the peak internal and external currents, 
plotted as function of Eo/Ecr- Field-current feedback be- 
comes important for Eo ^ 0.25 Ecr- 

system, so that the coupling between Eqs. (^) and (^) be- 
comes important. This marks the onset of field-current 
feedback, which can lead to plasma oscillations. 

Figure || depicts the ratio of the internal and external 
currents. Evidently, with the field strengths realisable at 
proposed XFELs, the internal current is negligible. How- 
ever, with the onset of particle accumulation the ratio 
grows with the number of elapsed laser periods and, since 
there are roughly 10^ periods in a given pulse, feedback 
will necessarily come to affect the plasma's behaviour at 
some stage in its evolution. Notably, for Eq = 0.5 Ecr 
that happens very early; i.e., after only 0.1 % of the life- 
time of the pulse, and feedback will have observable con- 
sequences. (NB. Our full calculations incorporate this ef- 
fect but the result is a signal that less complete estimates 
will become more unreliable in the domain of particle ac- 
cumulation.) A proper analysis of the observable effects 
of feedback must also include an accurate description of 



collisions, which can act to damp plasma oscillations [ [l5[ . 

We have explored the possibility of e~ e'^ pair produc- 
tion using XFEL beams as a parameter-free application 
of nonequilibrium quantum mean-field theory. An ide- 
alised model of the laser electric field was used as input 
to a quantum Vlasov equation coupled with Maxwell's 
equation. Decay of the QED vacuum is likely to be ob- 
served at proposed XFEL facilities. However, we showed 
that additional phenomena become observable if the peak 
power is increased by a factor of ~ 9. Then the particles 
produced from the QED vacuum accumulate to form a 
plasma, which is not achieved at the DESY and SLAG 
facilities. That plasma's behaviour will expose essential 
quantum field theoretical aspects of the pair production 
process, which are absent from the Schwinger formula. 
While this power is insufficient to expose negative en- 
ergy components in the electron wave packets, it can re- 
veal the non-Markovian nature of pair production; i.e., 
the feature that the appearance of a new pair is not in- 
dependent of the production of the preceding pairs nor 
the evolution of the plasma before this new pair is pro- 
duced. With this peak power the accumulated charged 
pairs are also likely to generate internal currents that 
interfere with the applied laser field so that the plasma 
exhibits charge-density fiuctuations; i.e., plasma oscilla- 
tions. An accurate prediction of this effect's observable 
consequences requires a realistic collision term in addi- 
tion to the proper source. 
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